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Abstract 

We investigate the prospects for the discovery of a neutral Higgs boson produced with one 
■ bottom quark fohowed by Higgs decay into a pair of bottom quarks at the CERN Large Hadron 

Colhder (LHC) and the Fermilab Tevatron ColHder. We work within the framework of the minimal 
supersymmetric standard model. The dominant physics background is calculated with realistic 
(~| I acceptance cuts and efficiencies including the production of bbb, bbb, jbb {j = g,q,q; q = u, d, s, c), 

tt — > bbjjiu, and tt — > bbjjjj. Promising results are found for the CP-odd pseudoscalar (A^) and 
CN ■ the heavier CP-even scalar (H^) Higgs bosons with masses up to 800 GeV for the LHC with an 

. integrated luminosity (L) of 30 fb^^ and up to 1 TeV for L = 300 fb~^. 
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I. INTRODUCTION 



The Fermilab Tevatron Run II has been taking data since March 2001, and the CERN 
Large Hadron Collider (LHC) is planned to start running in Autumn 2009. One of the 
most important experimental goals of the Tevatron Run II and the LHC is the search for 
the mechanism of electroweak symmetry breaking — to discover the Higgs bosons or to prove 
their non-existence. 

In the Standard Model, only one Higgs doublet is required to generate mass for both 
gauge bosons and elementary fermions, and the Higgs boson is the only particle remaining 
to be discovered in high energy experiments. In the minimal supersymmetric standard 
model (MSSM) [ij, the Higgs sector has Yukawa interactions with two doublets, 0i and 
02, whose neutral components couple to fermions with weak isospin = —1/2 and ts = 
+1/2 respectively 0]. After spontaneous symmetry breaking, there remain five physical 
Higgs bosons: a pair of singly charged Higgs bosons H^, two neutral CP-even scalars 
(heavier) and /i° (lighter), and a neutral CP-odd pseudoscalar A'^. The Higgs potential is 
constrained by supersymmetry such that all tree-level Higgs boson masses and couplings 
are determined by just two independent parameters, commonly chosen to be the mass of 
the CP-odd pseudoscalar (M^) and the ratio of vacuum expectation values of neutral Higgs 
fields (tan/3 = V2/V1). 

At the LHC, gluon fusion [gg — > 0°; 0° = h^,H^,A^) is the major source of neutral 
Higgs bosons in the MSSM for tan (3 less than 5. For tan^ > 7, neutral Higgs bosons are 
dominantly produced from bottom quark fusion bb 0° js], 0, [1, 0, 0] • Since the Yukawa 
couplings of (p^bb are enhanced by 1/ cos/5, the production rate of neutral Higgs bosons 
associated with bottom quarks, especially that of the or the H^, is enhanced at large 
tan p. 

For a Higgs boson produced along with a single bottorn quark at high transverse momen- 
tum (pt), the leading-order subprocess is bg — > 60° js], EO, 11, 12]. If two high px bottom 



quarks are required in association with a Higgs boson, the leading order subprocess should 



be gg bbcj) [3|, [131, [IJ, ll5|, Il6j. In 2002, it was suggested that the search at the LHC for 
a Higgs boson produced along with a single bottom quark with large px should be more 
pro mising than the production of a Higgs boson associated with two high px bottom quarks 
[lO|. This has already been shown to be the case for the decay mode of the Higgs 

bosons 17 . 



For large tan (3, the r r decay mode |18l . |19[ | can be a promising discovery channel for 



the A^ and the in the MSSM. Recently, the discovery channel bcfp 6r+r has been 
demonstrated to offer great promise at the LHC to search for the A^ and the H'^ up to 
Ma = 1 TeV [20j. 

The Higgs decay into bb has the largest branching fraction for large values of tan (3. 
However, the inclusive channel of pp —>■ (p^ ^ bb + X is very challenging at the LHC owing 
to the extremely large QCD background. Previous theoretical studies have focused on the 



associated production of bb(f) bbbb [2l|, |22, l23| . Realistic simulations by the ATLAS and 
the CMS collaborations with parton showering lead to pessimistic results jli, 25, 2^, because 
the trigger for the 46 final state requires high pt bottom quarks for pp — > 660" — > 6666 -|- X. 
The requirement of four high px 6-quarks removes most of the Higgs events. Moreover, 
integrating over the fourth 6-quark to study a 36 signal requires a careful inclusion of higher 
order corrections in the four-flavor scheme. These potentially large leading-log corrections 
are absorbed into the 6-quark PDFs in the five flavor scheme which we employ. 
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In this article, we present the prospects for discovering the MSSM neutral Higgs bosons 
produced with a single high pr bottom quark {b or b) followed by Higgs decay into a pair 
of bottom quarks . We calculate the Higgs signal and the dominant Standard Model (SM) 
backgrounds with exact matrix elements as well as realistic cuts and efficiencies. Further- 
more, we present promising 5a discovery contours at the LHC in the (Af4,tan/?) plane. 
Section II shows the production cross sections and branching fractions for the Higgs signal. 
The SM physics background is discussed in Section HI. Sections IV and V present the dis- 
covery potential at the LHC and the Fermilab Tevatron Run II. Optimistic conclusions are 
drawn in Section VI. 



II. THE PRODUCTION CROSS SECTIONS AND BRANCHING FRACTIONS 

At the LHC or the Tevatron Run II, the production cross section of bg — > b(f)'^ — »• bbb, 
where 0° = H^, h^, A^, is evaluated with the parton distribution functions of CTEQ6L1 [13] 
and the factorization scale iip = Mh/4: [l^. In this article, b represents a bottom quark 
(6) or a bottom anti-quark (6) unless it is explicitly specified. The bottom quark mass in 
the (p^bb Yukawa coupling is chosen to be the next-to-leading-order (NLO) running mass 
rrLhifiR) jlsl, which is calculated with mfe(pole) =4.7 GeV and the NLO evolution of the 
strong coupling ^]. We have also taken the renormalization scale to be Mh/4:. This choice 



of scale effectively reproduces the effects of next-to-leading order (NLO) corrections [10 
Therefore, we take the K factor to be one for the Higgs signal. 

At the LHC, we calculate the Higgs cross section cr(pp — ^ 60° bbb + X) with a Breit- 
Wigner resonance via bg — > 60° — > bbb. In addition, we check the cross section with the 
narrow width approximation 



cr 



{pp bcj)^ bbb + X) = a{pp bcj)^ + X) x 5(0° ^ bb) 



where 5(0° — > bb) is the branching fraction of a Higgs boson decay into bb. 

Figure 1 shows the invariant mass distribution (Mjj, i,j = 1,2,3) of the bibj or bibj pairs 
for the Higgs signal pp —>■ 6^4° — > bbb + X via bg — > 6^4°. The bottom quarks are ordered 
according to their transverse momenta, Pt(6i) > Pt(62) > Pribs). We note that with energy- 
momentum smearing, the cross section in the narrow width approximation (NWA) agrees 
very well with that evaluated via a Breit-Wigner resonance (BWR) for most parameters 
that we have chosen. Based on the ATLAS [24] specifications, we model these effects by 
Gaussian smearing of momenta: 

°-^°~0.03 (1) 



for jets at the LHC, with individual terms added in quadrature. For the Tevatron we use 

"■=''-0.03 (2) 



E y/E 

based on CDF parameters 3^. For M4 = 800 GeV and tan j3 = 50, the cross sections are in 
agreement within 10%. For large values of M^, the increased width of the Higgs may lead 
to a reduced signal due to cuts on the dijet invariant-mass acceptance window. This effect 
is less well-modeled in the NWA than with BWR, although the total cross-sections remain 
in good agreement. 
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FIG. 1: The invariant-mass distribution of 66 and 66 pairs da/dMbh{pp bbb + X), for the Higgs 
signal from bg 6A° with Ma = 200 GeV and tan /3 = 10 as well as Ma = 800 GeV for tan /3 = 10 
and tan /? = 50. We calculate the Higgs signal with a Breit-Wigner resonance (dash) and in the 
narrow width approximation (dot), applying minimal cuts oi px > 10 GeV and |?/| < 10. 



III. THE PHYSICS BACKGROUND 



The final state of bbb has dominant physics backgrounds coming from (a) bg — i 
(b) eg —>■ ebb, (c) qg — > qbb with q = u,d, s, (d) qq gbb with q = u,d, s, c, and (e) 
gg,qq —^tt^ bbjjiu, or gg,qq ^ tt ^ bbjjjj. We have computed the cross section of 
the Higgs signal and physics background utilizing MadGraph jsH, S^l and HELAS [ii to 



generate matrix elements. To reduce the physics background while keeping most of the 
signal events, we require that in each event there are three jets (at least two 6-jets) which 
satisfy the following requirements: 

(a) we consider two sets of cuts for an integrated luminosity (L) of 30 fb~^ (low luminosity, 
LL): (i) ptUi) > 50 GeV, pr(i2) > 30 GeV and ptOs) > 20 GeV (low pr cuts), or (ii) 
PtUi, j2, js) > 70 GeV (CMS 3-jet trigger) [1^ as well as the pseudorapidity, \r]\ < 2.5 
for all jets, where priji) > PT{j2) > PtUs), or 

(b) for L = 300 fb~^ (high luminosity, HL) we check two sets of cuts : (i) PriJi, j2, js) > 75 
GeV (ATLAS 3-jet trigger) Q or (ii) PtO'i, j2, js) > 150 GeV (ATLAS 3-jet trigger 
for high luminosity) [25| as well as \r]\ < 2.5 for all jets, 

(c) there is at least one pair of bottom quarks in the Higgs mass window such that 
\Mbb - M^l < AMhb, where AMbb = MAX(22 GeV, cta/), choosing ctm = 0.10 x or 
0.15 X for L = 30 fb"^ and cta/ = 0.15 x or 0.20 x for L = 300 fb'^, 

(d) all three jets are separated with AR = y/ Acp'^ + Ai]'^ > 0.7 (where is the angle 
between two jets in the transverse plane), 

(e) the missing transverse energy ( ^^) should be less than 40 GeV. 
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In addition, we veto events with more than three jets passing the cuts prij) > 15 GeV 
and \ri\ < 2.5. We take the 6— tagging efficiency to be eb = 0.6 (LL) or = 0.5 (HL), the 
probabihty that charm quark may be misidentified is ec = 0.15, and the probabihty that 
a hght quark or a gluon may be misidentified as a bottom quark is ej = 0.01. For the 
backgrounds arising from bbb and jbb [2l| as well as those from tt [sj], we assume a K factor 
of 2 when computing the significance as discussed below. In practice we find that the tt 
backgrounds are negligible after cuts, although we include them for completeness. 

In Figure 2, we present the transverse momentum distribution (da/dpT) of the bottom 
quarks {b or b), for the Higgs signal pp bbb + X. Also shown is the pt distribution 

for bottom quarks from the SM background bg bbb. We have required prib) > 10 GeV 
and |?7fe| < 2.5 in this figure. 

Vs = 14 TeV 
(a) bA" -> bbb", Mj = 200 GeV (b) bg -» bbb (SM) 
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FIG. 2: The transverse-momentum distribution for (a) the Higgs signal from bg — > bA^ with 
Ma_ = 200 GeV and tan /3 = 10, 50 as well as for (b) the physics background from bg — > bbb. We 
require prib) > 10 GeV and |%| < 2.5 in this figure. The vertical, dashed lines illustrate cuts at 
20 GeV and 70 GeV. 




IV. THE DISCOVERY POTENTIAL AT THE LHC 

To study the discovery potential of pp —>■ b(j)^ bbb + X [cjp = H^, h^, A^) at the LHC, 
we calculate the Higgs signal as well as the SM physics background in the mass window of 
M^±AMbb where AM^b = MAX(22 GeV, 0.10 x M^), or AMbb = MAX(22 GeV, 0.15 x M^) 
for an integrated luminosity of 30 fb~^. 

In Figure 3 we show the cross section of a{pp bbb + X), for tan (3 = 10 and 

50, with a common mass for scalar quarks, scalar leptons and the gluino fnj = '"^g = /i = 1 
TeV. We also present the background cross sections with no K factor in the mass window of 
Ma ± AMbb for the dominant SM processes pp bbb + X and pp jbb + X, j = q,q,g, with 
(a) low Pt cuts and (b) CMS 3-jet trigger. The cuts and tagging efficiencies are included 
with AMbb = 0.10 X M4. In addition, we present the 5a cross section for L = 30 fb~^. The 
cross section of the Higgs signal with tan/? ~ 50 can be larger than the 5a cross section 
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FIG. 3: The signal cross section of hg hA^ at the LHC for an integrated luminosity L = 30 
fb~^, as a function of Ma, for rriq = nig = fi = I TeV, tan /? = 10 and tan /? = 50. Also shown are 
the background cross sections in the mass window of Ma ±0.10 x Ma as discussed in the text for 
the SM contributions. We have applied acceptance cuts and efficiencies of tagging and mistagging. 



for Ma ^ 800 after acceptance cuts. Requiring higher transverse momenta {px > 70 GeV) 
greatly reduces the background and the Higgs signal for Ma < 200 GeV. 

We define the signal to be observable if the lower limit on the signal plus background is 



larger than the corresponding upper limit on the background [35|, |36j, namely. 



L{as + (Tb) - NJL{as + (^b) > Lab + NJLa, 



which corresponds to 



a,, > 



N_ 

T 



2 r 



1 + 2JLab/N 



(3) 



(4) 



Here L is the integrated luminosity, as is the cross section of the Higgs signal, and ab is 
the background cross section. Both cross sections are taken to be within a bin of width 
±AMbb centered at M^. In this convention, N = 2.5 corresponds to a 5a signal. We take 
the integrated luminosity L to be 30 fb~^ and 300 fb~^ (23 |. 

For tan/3 > 10, Ma and Mh are almost degenerate when Ma ^ 125 GeV, while Ma and 
m/j are very close to each other for Ma ^ 125 GeV in the MSSM Therefore, when 

computing the discovery reach, we add the cross sections of the and the /i" for Ma < 125 



38 



GeV and those of the A° and the for Ma > 125 GeV [2J, [2 

Figure 4 shows the 5a discovery contours for the MSSM Higgs bosons where the discovery 
region is the part of the parameter space above the contour. We have chosen Msusy = ''^q = 
nig = mi = ^ = 1 TeV. If Msusy is smaller, the discovery region of A^, — > bb will be 
slightly reduced for Ma ^ 250 GeV, because the Higgs bosons can decay into supersymmetric 
(SUSY) particles [3^ and the branching fraction of 0° —* bb is suppressed. For Ma ^ 125 
GeV, the discovery region of bb is slightly enlarged for a smaller Msusy? but the 

observable region of —>■ bb is slightly reduced because the lighter top squarks make the 
and the lighter; also the H'^bb coupling is enhanced while the h'^bb coupling is reduced (38 
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FIG. 4: The 5a discovery contours at the LHC with ^/s = 14 TeV for (a) L = 30 fb ^ and low px 
cuts, (b) L = 30 fb-i and pr > 70 GeV, (c) L = 300 fb^^ and pr > 75 GeV, (d) L = 300 fb"! and 
Pt > 150 GeV, in the Ma versus tan /3 plane. The signal includes (j)^ = and for Ma < 125 
GeV, and (/>° = and for Ma > 125 GeV. The discovery region is the part of the parameter 
space above the contours. 

In addition, we have studied the effect of an invariant mass cut, using only the two jets 
with highest pr as the candidate pair. Table I presents the cross section corresponding to 
two schemes: (a) requiring \Mi2 — M^\ < AMf,f,, and (b) requiring |Mjj — M^l < AMf,f,;i,j = 
1, 2, 3. We find that for Ma ^ 400, it is more advantageous to apply an invariant mass cut 
only on the two leading b jets. For lower masses using any pair of the three leading jets 
leads to higher significance. We also show the ratio of signal to background in this figure. 
We have chosen a set of cuts, PtUi, j2, js) > 100, 80, 70 GeV, which tends to maximize this 
ratio. Less stringent cuts can improve the nominal statistical significance in the low mass 
regions, but for high masses and low tan f3 the small signal to background ratio would require 
excellent understanding of backgrounds and systematic errors. 

Furthermore, we have studied the effects of SUSY particles on the (p^bb Yukawa couplings 
at large tan The SUSY contributions can be described with an effective Lagrangian and a 



7 



TABLE I: Cross sections in fb for the Higgs signal and physics background for two choices of 
cuts on the invariant mass of bb: (a) two leading jets (M12) versus (b) any two jets {Mij) used to 
reconstruct the Higgs invariant mass. Significances are computed with L = 30 fb^^. 



Mass(GeV) 




Signal 


Background 


Significance {Nss = Ns/V^b + Ns) 


Ns/Nb 


tan P =10 


200 


M12 
M,, 


44.2 
126 


3960 
14500 


3.82 
5.72 


1.12 X 10~2 
8.70 X 10-3 


400 


M12 
Mij 


23.5 
32.2 


6680 
11900 


1.57 
1.61 


3.52 X 10-3 
2.70 X 10-3 


800 


M12 
Mij 


1.42 
1.61 


1400 
2380 


0.208 
0.181 


1.02 X 10-3 
6.76 X 10"^ 


tan P = 20 


200 


M12 
Mi, 


178 
498 


3960 
14500 


15.1 
22.2 


4.48 X 10-2 
3.43 X 10-2 


400 


M12 
Mi, 


104 
143 


6680 
11900 


6.94 
7.14 


1.56 X 10-2 
1.20 X 10-2 


800 


M12 
Mi, 


6.99 
7.96 


1400 
2380 


1.02 
0.891 


5.00 X 10-3 
3.34 X 10-3 


tan/? = 50 


200 


A/12 
Mi, 


961 
2770 


3960 
14500 


75.0 
115 


2.42 X 10-1 
1.91 X 10-1 


400 


M12 
Mi, 


563 
792 


6680 
11900 


36.2 
38.5 


8.43 X 10-2 
6.66 X 10-2 


800 


M12 
Mi, 


38.7 
44.7 


1400 
2380 


5.58 
4.96 


2.76 X 10-2 
1.87 X 10-2 



function Ab 4^, 41, H, 43 1 such that the bottom quark mass in Yukawa couphngs becomes 

nib 



nib 



where SUSY QCD corrections lead to 

^Jm-giJ, tan /31{mi^ , m^^ , m-g) 



for bottom squarks and gluinos, and the auxihary function is 

1 



/(a, b, c) 



a%Hn^ + b'c'\n^- + c'a'\n-) 



6_2 

c2 



(a2-62)(52_c2)(c2-a2) 
Then the cross section can be estimated with a simple formula 

tan2/5 



c 

^2 



a{pp h(f + X)xB{ 



bb) - asAiipp ^ bH + X) X 
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(1 + Ab) 



X 



:i + Afe)2 + 9 



TABLE II: Effect of Af, in M^"^ (no mixing) scenario. Cross sections in fb for pp — > b<p^ bbb+X 
using high px (> 70 GeV) cuts. Tagging efficiencies have not been applied. 



Mass(GeV) 


Ab = 


g/b 


g/b + H/t 


tan P = 10 


200 


= 


+200 


698(708) 


646(658) 


633(656) 




= 


-200 


699(703) 


745(755) 


761(753) 


400 


= 


+200 


155(155) 


143(144) 


140(145) 




= 


-200 


156(155) 


168(167) 


172(168) 


800 


= 


+200 


7.90(7.91) 


7.28(7.31) 


7.07(7.31) 




= 


-200 


7.87(7.93) 


8.63(8.56) 


8.86(8.60) 


tan/3 = 50 


200 


= 


+200 


16400(16400) 


12200(12200) 


11000(12200) 




A* = 


-200 


16400(16300) 


22600(22600) 


25800(22600) 


400 


= 


+200 


4120(4120) 


3060(3060) 


2750(3060) 




At = 


-200 


4130(4120) 


5730(5730) 


6560(5720) 


800 


= 


+200 


233(233) 


172(172) 


154(172) 




= 


-200 


233(233) 


325(325) 


373(325) 



In our analysis of SUSY effects, we adopt the conventions in Refs. |12l . |4J] and have used a 
more complete estimate, including the effects of the modified Higgs width in the full BWR 
calculation. Table II shows the cross section including (a) no SUSY effects, (b) contributions 
from bottom squarks and gluinos, and (c) contributions from bottom squarks and gluinos as 
well as from top squarks and Higgsinos. The top squark/Higgsino loops give an additional 
effective correction to rrib, 



Afe = — At fi tan (31 {mt^,m^^, fi) 



where at = A^/47r (At = y/2'mt/v2 being the top Yukawa coupling), and At is the trilinear 
Higgs-stop coupling. It is clear that SUSY effects reduce the Higgs cross section for a positive 
/i while they enhance the Higgs cross section for a negative fi. The effect of the Higgsino/stop 
loops is highly dependant on the size of At. We present two scenarios, M^"^ and no-mixing, 
as defined in Ref. In the former the Higgsino/stop contribution is comparable to the 
gluino/bottom-squark term, in the latter it is almost negligible. 



V. THE DISCOVERY POTENTIAL AT THE FERMILAB TEVATRON 



To study the discovery potential of Higgs decays into bottom quark pairs at the Fermilab 
Tevatron Run II, we require 

(i) three b quarks or 3 jets (at least two b jets) with px > 15 GeV or PxiJi, j2, js) > 
50,30, 15 GeV, \ri{b,j)\ < 2.0, and a 6— tagging efficiency = 50% |30l ]. 



(ii) the angular separation between each pair of jets should be AR > 0.4 [45|, 
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(iii) the invariant mass of the reconstructed bottom quark pairs should be within the mass 
window of the Higgs mass with AMbb = MAX(0.1 x M^, 20GeV). 

Figure 5 show the 5a discovery contours for the MSSM Higgs bosons, where the discovery 
region is the part of the parameter space above the curves. The discovery contours for 
AMbb = 0.10 X Iii are comparable to those presented in this figure. 

We find that the Tevatron Run II can discover neutral Higgs bosons in the MSSM for a 
value of tan/3 slightly less than 30 with an integrated luminosity of 4 fb~^ and Ma < 120 
GeV. For tan /? ~ 50, the Tevatron Run II will be able to discovery the Higgs bosons up to 
Ma ~ 160 GeV with L = 4 fb-\ and up to _M^_~ 200 GeV with L = 20 fb'^ Our results 
are consistent with those found in Refs. [ 



45, 47 
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FIG. 5: The 5(t discovery contours at the Fermilab Tevatron Run II for an integrated luminosity 
(L) of 4 fb~\ 10 fb~\ 20 fb~^ in the Ma versus tan/3 plane. The signal includes = and /i° 
for Ma < 125 GeV, and (p^ = and for Ma > 125 GeV. The discovery region is the part of 
the parameter space above the contours. 



VI. CONCLUSIONS 

The associated production of a Higgs boson with a bottom quark, followed by the Higgs 
decay into bottom quark pairs, is a promising channel for the discovery of the neutral Higgs 
bosons in the minimal supersymmetric standard model at the LHC. The and the 
should be observable in a large region of parameter space with tan/3 ^ 10. The associated 
final state of 60° bbb could discover the and the at the LHC with an integrated 
luminosity of 30 fb~^ if Ma ~ 800 GeV. At a higher luminosity of 300 fb^^, the discovery 
region in Ma is expanded up to Ma = 1 TeV for tan/3 ~ 50. 

In our analysis, we apply a mass cut, requiring the reconstructed Higgs mass to lie in 
the mass window M^ ± AMbb- We note that improvements in the discovery potential will 
be possible by narrowing AMbb if the bottom quark pair mass resolution can be improved. 
In regions of high mass and low tan/3 the ratio of signal to background events is very 
low. Discovery in these regions would require either excellent understanding of backgrounds 
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in order to lower systematic errors below the few percent level, or better discrimination 
between signal and background due to narrower AMbb or improved 6-tagging. Our results 
using three 6's are more promising than those found in previous studies based on 46 analyses 
2l|,[22l,[2l,[23. 



The discovery of the associated final state of 60° 
b(f) b^^fi' 



bbb along with 60° br^r 



[20| and 

17| will provide information about the Yukawa couplings of //0°; f = b,T, fj,, 
for fermions with ^3 = —1/2. Furthermore, the muon pair channel can also be observable in 
a significantly large region and the muon pair channel will pro vide a good opportunity to 
precisely reconstruct the masses for MSSM Higgs bosons [l^,]l7, 3^. In concert, this family 
of channels may provide an excellent window on the Yukawa sector of the MSSM. 
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